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In  this  paper  we  consider  the  complete  thermofluid  design  and  performance  of  a  thermoelectric  module. 
We  increase  the  temperature  difference  that  must  be  maintained  across  the  module,  and  at  the  same 
time  we  reduce  the  pumping  power  required  by  the  streams  that  bathe  the  hot  and  cold  plates  of  the 
module.  We  find  that  for  greater  power  output  the  two  streams  must  be  configured  in  parallel,  not  in 
counterflow,  and  not  between  two  well  mixed  plenums.  We  also  find  that  the  loss  of  thermoelectric 
power  due  to  the  temperature  nonuniformity  of  the  two  plates  competes  with  the  power  lost  during 
the  pumping  of  the  two  streams,  and  that  from  this  competition  results  the  optimal  mass  flow  rates  of 
the  two  streams.  At  the  optimum,  the  maximized  power  output  of  the  module  is  proportional  to  a  group 
of  geometric  parameters  (Eq.  (39)),  which  can  be  maximized  further  by  designing  vascular  flow  architec¬ 
tures  for  the  two  plates.  The  vascular  designs  reveal  an  optimal  ratio  of  channel  diameters,  optimal  plate 
aspect  ratio,  and  a  channel  flow  volume  fraction  that  is  of  the  same  order  as  the  ratio  of  the  fluid  thermal 
conductivity  divided  by  the  solid  thermal  conductivity.  The  flow  architectures  are  further  illustrated  with 
numerical  examples  of  vascular  and  serpentine  configuration  and  performance. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Thermoelectric  power  generation 

Thermoelectric  power  generation  (TE)  is  attracting  interest 
from  across  the  spectrum  of  technology  because  it  offers  a  path 
to  clean,  renewable  energy.  In  addition,  TE  devices  offer  compact¬ 
ness  and  simplicity  because  they  do  not  have  moving  parts.  The 
efficiency  of  TE  designs  is  limited  by  the  transport  properties  of 
the  thermoelectric  materials  accounted  for  in  the  property  group 
Z  discussed  in  Eq.  ( 1 ).  Progress  toward  greater  efficiencies  can  be 
made  by  optimizing  the  geometry— the  relative  shapes  and 
sizes— of  the  legs  of  the  module  [1-3].  This  and  other  aspects  of 
the  optimization  of  thermoelectric  modules  (TEM)  has  generated 
a  growing  body  of  literature  [1-19]. 

In  the  present  study  we  consider  the  whole  picture— the  ther¬ 
mofluid  design  and  performance  of  the  TEM.  We  consider  not  only 
the  temperature  difference  that  must  be  maintained  across  the 
TEM,  but  also  the  fluid  mechanics  of  the  vascular  hot  and  cold 
plates  that  sandwich  the  TEM.  We  show  that  important  tradeoffs 
exist  in  the  thermofluid  design  such  that  the  overall  performance 
of  the  TEM  is  increased. 

The  thermoelectric  module  (TEM)  is  defined  in  Fig.  1.  The  hot 
plate  is  heated  by  a  stream  of  hot  oil.  The  cold  plate  is  cooled  by 
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a  stream  of  water.  Heat  flows  by  conduction  perpendicular  to  the 
plates,  from  the  hot  plate  to  the  cold  plate. 

Larger  stacks  are  obtained  by  sandwiching  several  TEMs  as 
shown  in  Fig.  2.  All  the  TEMs  are  identical.  Each  is  a  sandwich  of 
three  parts:  hot  plate  (H),  thermoelectric  converter  (T),  and  cold 
plate  (C),  in  other  words,  HTC,  or  CTH.  There  are  two  ways  to 
assemble  such  elements  in  a  stack: 

(a)  Alternating  orientations,  i.e.,  plates  of  the  same  temperature 
touching  each  other: 

HTC  CTH  HTC... 

(b)  Elements  oriented  in  the  same  way: 

HTC  I  HTC  I  HTC... 

The  simplest  design  is  (a).  In  design  (b),  a  layer  of  insulation  (I)  must 
be  placed  between  the  C  and  H  plates  of  adjacent  elements.  Fig.  2 
shows  the  (a)  rule  of  assembly. 

If  the  TEM  is  small  enough  and  the  hot  and  cold  fluid  flow  rates 
are  large  enough,  then  the  hot  and  cold  plates  are  essentially  iso¬ 
thermal,  at  Th  and  TL,  respectively.  The  thermodynamics  of  the 
TEM  is  detailed  in  Ref.  [1  ],  and  it  shows  that  the  maximum  thermo¬ 
dynamic  efficiency  of  a  module  with  isothermal  ends  is 
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Nomenclature 

A 

plate  area,  m2 

W 

power,  W 

C 

factor,  dimensionless 

Wiost 

lost  power,  W 

Cp 

specific  heat  at  constant  pressure,  J  kg  1  K  1 

Wmax 

maximum  power,  W 

d 

spacing  between  channels,  m 

WAP 

pumping  power,  W 

d2 

channel  thickness  scale,  m 

WAT 

power  lost  due  to  temperature  nonuniformity,  W 

Dr 

trunk  thickness,  m 

X 

longitudinal  coordinate,  m 

d2 

branch  thickness,  m 

z 

group  of  properties,  Eq.  (1) 

G 

function  of  geometry,  Eq.  (39) 

z 

function  of  Z,  Eq.  (26) 

H 

transversal  dimension  of  area,  m 

I 

integral 

Greek  symbols 

k 

effective  TEM  thermal  conductivity,  W  m_1  K  1 

S 

spacing  between  branches,  m 

kf 

fluid  thermal  conductivity,  W  m_1  K  1 

A  P 

pressure  difference,  Pa 

ks 

solid  plate  conductivity,  W  m  1  K  1 

AT 

temperature  difference,  K 

Kl.2 

factors,  Eqs.  (37)  and  (38) 

n 

efficiency,  Eq.  (2) 

L 

transverse  dimension  of  area,  m 

RP 

pump  isentropic  efficiency 

Ls 

length  scale,  m 

kinematic  viscosity,  m2  s  1 

m 

number  of  channels  (branches) 

{ 

dimensionless  longitudinal  coordinate 

m 

mass  flow  rate,  kg  s  1 

p 

density,  kg  m  3 

n 

number  of  channels  (branches) 

<t> 

area  fraction 

N 

number  of  heat  transfer  units,  Eq.  (1 1 ) 

<p 

porosity 

P 

perimeter  of  contact,  m 

<t>c 

critical  area  fraction 

P 

pressure,  Pa 

Q  H 

R 

heat  current,  W 
electric  resistance,  Cl 

Subscripts 

Rth 

thermal  resistance,  K  W  1 

fluid 

Sv 

t 

svelteness  number 
thickness  of  TEM,  m 

P 

counter 

counter  flow 

t 

spacing  between  channels,  m 

h 

hot 

Th 

high  temperature,  K 

in 

Tl 

low  temperature,  K 

minimum 

TEM 

thermoelectric  module 

opt 

para 

optimum 
parallel  flow 

U 

V 

overall  heat  transfer  coefficient,  W  m  2  K  ' 
volume,  m3 

Th-Tl  Z-l 
Th  Z+Tl/Th 


(1) 


where  Z  (constant)  represents  a  group  of  thermophysical  properties 
of  the  material  structure  sandwiched  between  TH  and  TL.  The  energy 
conversion  efficiency  is  defined  as 


n 


Qh 


(2) 


where  I2R  is  the  delivered  electrical  power,  and  qH  is  the  total  heat 
transfer  rate  from  TH  to  TL.  The  power  I2R  is  associated  with  the  flow 
of  the  generated  current  I  through  an  external  electrical  resistance 
R. 

The  heat  current  qH  is  proportional  to  the  temperature  differ¬ 
ence  that  drives  it, 


<7h§. 


H:-  Tl 
5th 


(3) 


where  Rth  is  the  thermal  resistance  of  the  TEM  structure  held 
between  TH  and  TL.  By  combining  Eqs.  (1 )— (3)  we  arrive  at  the 
important  conclusion  that  maximum  power  is  delivered  in  propor¬ 
tion  with  (Th  —  TL)2, 

(f2R)max  =  (Th  -  Tl)2  •  Rth7-H(^1TL/rH)  (4) 

The  chief  conclusion  that  follows  from  thermodynamics  [1]  is 
that  maximum  power  is  achieved  in  proportion  with  the  tempera¬ 
ture  difference  squared.  When  the  volume  of  the  TEM  is  fixed, 
maximum  power  also  means  maximum  power  density.  The  route 


to  better  designs  with  greater  power  density  consists  of  maintain¬ 
ing  the  largest  temperature  difference  across  the  TEM,  between 
any  point  on  the  hot  plate  and  its  projection  on  the  cold  plate. 

When  the  hot  and  cold  plates  are  not  isothermal,  the  design  rule 
identified  above  becomes  critical,  and  it  applies  at  every  point  in 
the  plane  of  the  hot  plate  or  the  cold  plate. 

The  largest  temperature  difference  available  is  the  difference 
between  the  inlet  of  the  hot  stream  (Th,m)  and  the  inlet  of  the  cold 
stream  (Tc  in).  Ideally,  the  thermo-fluid  design  of  the  TEM  should  be 
such  that  the  hot  plate  is  isothermal  at  rh,in,  and  the  cold  plate  is 
isothermal  at  Tc  in.  This  design  would  require  thick  high-conductiv- 
ity  bodies  for  the  hot  and  cold  plates,  which  increase  the  weight 
and  decrease  the  power  density  of  the  TEM. 

The  challenge  is  to  design  the  TEM  when  its  hot  and  cold  plates 
are  large  and  not  isothermal.  This  can  be  achieved  by  vascularizing 
the  two  plates  such  that,  although  nonuniform,  the  local  tempera¬ 
ture  difference  (Th  -  Tc)  distribution  is  such  that  the  integral  of 
(Th  -  Tc)2  over  the  plane  of  the  plate  is  the  largest  (Fig.  3). 


2.  Parallel  flow,  or  counterflow? 

Consider  first  the  unidirectional  flow  configuration  shown  in 
Fig.  4.  The  hot  stream  flows  as  a  sheet  in  one  direction  (x)  parallel 
to  the  sheet  flow  of  cold  fluid.  The  capacity  rates  of  the  two  streams 
are  the  same,  mcP.  The  stream  temperatures  vary  in  the  longitudi¬ 
nal  direction,  Th(x)  and  Tc(x).  When  the  streams  are  configured  in 
parallel  flow  (Fig.  4),  the  first  law  of  thermodynamics  for  each 
stream  control  volume  of  length  dx  requires 
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Fig.  1.  The  main  features  of  the  thermoelectric  module. 


-mcPdTh  =  q'dx  (5) 

mcPdTc  =  q'dx  (6) 

where  the  stream-to-stream  heat  transfer  rate  q’  [W/m]  is 

q'  =  Up(Th  -  Tc)  (7) 

where  p  is  the  line  of  thermal  contact  in  the  direction  perpendicular 
to  Fig.  4,  and  U  =  k/t  is  a  constant  representing  the  effective  thermal 
conductivity  (k)  between  the  Th  and  Tc  plates,  and  the  thickness  of 
the  structure  held  between  the  two  plates  (t).  Combining  Eqs. 
(5)— (7)  and  integrating  from  the  common  inlet  (x  =  0)  where  Th- 
—  Tc=  Thin  -  Tcin,  we  obtain  the  temperature  difference  distribution 


AT  =  ATine~m 

(8) 

where  AT  =  Th(x)  -  Tc(x) 

(9) 

ATj„  =  Th  i„  -  TC  i„ 

(10) 

N  =  ^ 

(11) 

mcp 

H 

(12) 

and  pL  =  A  is  the  area  of  the  plate,  in  other  words,  N=  (k/t)A/(mCp). 
According  to  Eq.  (4),  the  power  density  of  the  parallel  flow  design 
is  proportional  to  the  integral 
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Tc,i„  ■ 


I _ >  rhcp  Tc(x) 

0  x 

Fig.  4.  Two  parallel  streams  flowing  in  the  san 


W  =  /V2dc  =  ^(1- e-2' 


(13) 


available  work)  in  heat  exchangers  demonstrates  that  counterflow 
is  a  better  design  than  parallel  flow  (e.g.  Ref.  [1,  p.  602]). 

Why  both  conclusions  are  correct,  parallel  flow  for  the  present 
TEM  objective,  and  counterflow  for  power  plant  efficiency  broadly 
Ref.  [1,  p.  602]  is  thermodynamics.  Key  in  the  present  paper  is  that 
the  parallel  flow  offers  greater  and  greater  power  densities  relative 
to  counterflow  as  N  increases,  Eq.  (19),  that  is  as  the  fluid  spends 
more  time  in  residence  in  its  channels,  in  longer  and  more  tortuous 
channels.  In  this  friction  dominated  limit,  in  addition  to  pursuing 
greater  power  density  of  generation  it  is  important  to  pursue  less 
and  less  pumping  power  for  the  fluid  streams  that  bathe  the 
two-flow  structure  of  Fig.  4. 


When  the  two  streams  are  arranged  in  counterflow,  the  local 
stream-to-stream  temperature  difference  is  independent  of  x.  The 
counterflow  analysis  that  replaces  Eqs.  (5)  and  (6)  begins  with 

-mcpdTh  =  q'dx  (14) 

-  mcp  dTc  =  q'dx  (15) 


3.  Mixed  plenums 

Another  design  alternative  is  to  fit  the  hot  end  and  cold  end  of 
the  TEM  with  plenums  in  which  the  fluid  is  well  mixed  at  TH  and  TL, 
respectively.  In  the  hot  plenum,  the  fluid  temperature  drops  from 
Thjn  to  Th,  while  in  the  cold  plenum  it  rises  from  Tcin  to  TL, 


and,  after  using  Eq.  (7),  shows  that  Eq.  (8)  is  replaced  by 

AT-  ATi" 

AT  1 +N 


The  power  density  integral  (13)  is  replaced  by 


a  r?„ 
(1+JVf 


(16) 


(17) 


mcP  (Th,in  -  Th)  =  UpL(T„  -  TL) 
mcp(TL-TcM)  =  UpL(TH-TL) 

Adding  Eqs.  (20)  and  (21),  we  find  that 
AT  A  Tin 

AT  1  +  2N 

with  the  corresponding  performance  integral 


(20) 

(21) 


(22) 


The  relative  merit  of  the  parallel  flow  and  the  counterflow  con¬ 
figurations  is  indicated  by  the  ratio  of  Eqs.  (13)  and  (17): 


(23) 


(1  +  N)2 

2N 


(1  -  e~2N) 


(18) 


When  N  »  1,  this  ratio  approaches  a  value  much  greater  than  1, 


(19) 


When  N«l,  Eq.  (18)  approaches  1.  In  between,  for  example  at 
N=  1,  Eq.  (18)  approaches  1.264,  which  indicates  that  the  power 
density  of  the  parallel  flow  configuration  exceeds  by  26.4%  the 
power  density  of  the  counterflow  configuration. 

The  conclusion  that  parallel  flow  is  more  advantageous  than 
counterflow  is  important  and  fundamental.  The  fine  channels  with 
parallel  flows  are  shown  in  Fig.  5.  This  finding  is  important  because 
contrary  to  the  first-designs  shown  in  Figs.  1  -3,  the  flow  in  the  fine 
transversal  flow  channels  should  be  oriented  in  parallel  flow,  not 
counterflow.  It  is  fundamental  because  all  the  literature  on  the 
minimization  of  irreversibility  (entropy  generation,  loss  of 


The  above  integral  is  smaller  than  the  integral  for  counterflow, 


/counter  _  (1  +-3M f 
/mixed  ~~  (1  +N)2 


(24) 


This  ratio  tends  to  4  as  N  »  1.  In  summary,  the  relative  performance 
of  the  three  designs  is  in  favor  of  the  parallel  flow  configuration,  cf. 
Fig.  6, 

/para  >  /counter  >  /mixed  (25) 

and  the  strength  of  the  inequality  signs  increases  monotonically  as 
N  increases.  In  the  opposite  limit,  N<  1,  the  inequality  signs  are 
replaced  by  equal  signs. 

Friction  is  also  a  detrimental  feature  in  the  limit  of  fast  flow 
where  N  -c  1  and  the  two  plates  approach  the  isothermal  descrip¬ 
tion.  This  further  enhances  the  importance  of  designing  the  flow 
architecture  of  each  plate.  We  focus  on  this  feature  of  the  design 
in  the  next  two  sections. 
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N 

Fig.  6.  The  power  performance  of  parallel  flow,  counter  flow  and  mixed  plenums. 


4.  Heat  transfer  and  fluid  friction  losses 

Consider  next  the  broad  view  of  the  thermodynamic  perfor¬ 
mance  of  the  TEM.  The  uppermost  level  of  performance  is  indi¬ 
cated  in  Eq.  (4),  which  can  be  rewritten  as 

Wjnax  =  ATfnZ  (26) 

where  Wmax  is  the  maximum  power  output  (/2P)max,  ATin  =  TH  -  TL, 
and  Z  is  the  Z-dependent  function  appearing  on  the  right  side  of  Eq. 
(4).  The  actual  power  output  (W)  is  smaller  than  Wmax  because  of 
two  unavoidable  design  features:  the  temperature  distribution  on 
the  two  plates  is  not  perfectly  isothermal,  and  pumping  power  must 
be  spent  in  order  to  push  the  hot  and  cold  fluids  through  the  vascu¬ 
latures  built  into  the  plates. 

Assume  that  the  orientation  of  the  two  streams  is  in  parallel 
flow,  because  according  to  the  inequalities  (25)  the  parallel  flow 
design  is  superior  to  the  other  designs  for  any  finite  N.  In  the 
absence  of  pumping  power  losses,  the  actual  power  output  W  is 
equal  to  fpara  Z,  cf.  Eqs.  (4)  and  (13).  When  the  pumping  power 
Wap  is  not  negligible,  the  actual  power  output  is 


which  becomes 

Wlost  =  WAT  +  W^  (29) 

where  the  WAr  term  accounts  for  the  loss  due  to  the  departure  from 
temperature  uniformity  on  the  two  plates, 

WAP  =  ATfnz(l-l^N)  (30) 

The  dependence  of  WAT  on  N  is  illustrated  in  Fig.  7 :  the  heat  transfer 
loss  decreases  monotonically  as  N  decreases,  and  in  the  limit  JV  <g  1 
it  is  represented  by 

Ww^ATlZN  (31) 

The  limit  N  c  1  means  that  the  capacity  rate  mcP  is  high  enough  to 
overtake  (k/t)A  in  the  N  definition  (11),  namely  N  =  (|)A/(mcP)  In 
this  limit  the  flow  rate  m  may  be  high  enough  so  that  the  pumping 
power  loss  Wap  is  not  negligible  relative  to  WAr. 

The  pump  power  is  Wap  =  2mAP/(pi]p)  where  >/p  ^  1  is  the 
pumping  isentropic  efficiency,  assumed  constant.  The  factor  2 
means  that  there  are  two  streams  that  must  be  pumped,  i.e.  two 
plates  with  channels. 

Next,  assume  that  the  flow  is  fully  developed  and  laminar  in  all 
the  channels  of  the  plate  vasculature,  and  that  the  local  AP  losses 
are  negligible.  In  this  regime,  the  overall  AP  that  must  be  main¬ 
tained  between  the  inlet  and  outlet  of  one  plate  obeys  the  form 
[20] 

A P  =  Cvm±  (32) 

where  Ls  and  D*  are  the  scales  (dimensions)  of  the  flow  path  length 
and  channel  diameters,  and  C  is  a  dimensionless  constant.  In  the 
next  section  we  show  that  Ls,  Ds  and  C  depend  on  the  configuration 
and  size  of  the  plate  vasculature.  At  this  point,  we  retain  the  conclu¬ 
sion  that  the  overall  pumping  power  increases  monotonically  with 


Wap  =  2—  ^m2  (33) 

P%D4S 

Adding  WAp  to  the  heat  transfer  loss  in  the  small-N  limit,  Eq.  (31), 
we  obtain  the  total  power  loss  expression 


W  =  fparaZ  -  Wap  (27) 

The  lost  power  (W|ost  =  Wmax  -  W)  is  due  to  heat  transfer  and  fluid 
friction  effects  combined, 


Wlost  =  AT2,Z-IparaZ+WAP 


(28) 


n  Cv  Ls  .  t 

+  2 - =^m2 

P%D*S 


(34) 


which  is  of  the  form  Wlost  =  K,rh  1  +  K2m2.  Fig.  8  shows  that  Wk 
has  a  minimum,  which  is  located  at 


N 


Fig.  7.  The  power  loss  W&T  decreases 


lically  as  N  decreases. 
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Fig.  8.  The  mass  flow  rate  m  resulting  from  the  tradeoff  between  the  AT  and  AP 


Wlosmin  =  ^K]'3Kl'3 
with  the  notation 


2  Cv  Is 

WpDf 


(35) 

(36) 


(37) 

(38) 


Noteworthy  is  the  effect  of  TEM  geometry:  Kj  is  proportional  to  A(t, 
and  K2  is  proportional  to  CLJD4.  This  means  that  the  minimum  total 
loss  of  power  is  proportional  to  the  geometric  factor 


Further  reductions  in  WW  min  are  possible  by  seeking  flow  architec¬ 
tures  with  smaller  G  values  subject  to  the  volumetric  constraints 
that  govern  the  TEM.  This  aspect  of  the  constructal  design  explora¬ 
tion  is  detailed  in  the  next  section. 


5.  Vascular  architecture 


A  P1  =  Cvifh^  (40) 

where  the  factor  i  is  due  to  the  linear  variation  of  the  mass  flow 
rate,  and  C  is  the  constant  noted  in  Eq.  (32),  for  example,  C  =  128/ 
n  if  the  channel  cross  section  is  round. 

Next,  assume  that  the  m  stream  is  distributed  almost  uniformly 
to  the  n  transversal  branches.  The  applicability  of  this  approxima¬ 
tion  is  treated  in  detail  in  Ref.  [22],  The  pressure  drop  along  one  D2 
channel  is 

AP2  =  (41) 

The  overall  pressure  difference  that  drives  the  rh  stream  is 
AP  =  AP,  +  A P2.  The  objective  of  the  following  analysis  is  to  search 
for  the  flow  configuration  with  the  smallest  flow  resistance  (A P/m) 
subject  to  the  global  finite-size  constraints  of  the  design.  Combining 
Eqs.  (40)  and  (41 )  leads  to 
AP  1  ,  H 

Cvm  2D\ +  nD42  (  ) 

which  can  vary  subject  to  A  =  HL.  The  flow  resistance  is  minimum 


H  =  n/DA4 

1  2  VdJ 

(43) 

r  (2 A\ 1/2 /DA  2 

(44) 

W  (d,) 

(45) 

AP™,  (2 A/n)1/2 

Cvm  ~  (D,D2)2 

(46) 

Further  reductions  in  flow  resistance  can  be  achieved  by  maximiz¬ 
ing  the  product  D,D2.  The  two  dimensions  (D,,D2)  are  related 
through  the  requirement  that  the  “printed  area”  of  the  vasculature 
must  be  a  fraction  (</>)  of  the  plate  area: 

</>  =  i  (2D,L  +  nD2H)  (47) 


Consider  the  flow  architecture  embedded  in  one  plate,  Fig.  9, 
which  consists  of  two  trees  matched  canopy  to  canopy  [20],  The 
first  tree  distributes  the  m  stream  to  the  plate  area  A  =  HL,  and 
the  second  collects  the  stream  from  the  area  and  guides  it  to  the 
outlet.  The  trunk  is  a  channel  with  transversal  dimension  D,,  and 
the  n  branches  have  the  dimension  D2.  Assume  that  n  is  large  such 
that  along  the  trunk  the  flow  rate  varies  almost  linearly  between 
the  opening  (inlet,  or  outlet)  and  the  closed  end.  In  this  case,  the 
pressure  drop  along  one  of  the  D,  trunks  is 


Fig.  9.  Vascular  design  consisting  of  two  trees  matched  canopy  to  canopy. 


By  using  Eqs.  (44)  and  (45),  the  </>  constraint  becomes 
(2A),/2_4f^ 

^  t**®  =n2D\  +  D\ 

Minimizing  the  expression  (46)  subject  to  the  constraint 
obtain 


(48) 
(48)  we 


°I  =  mv* 

D2  \2 ) 

(49) 

=  </>rr3/102-17/10 

A' 

(50) 

=  </m-7/102-,3/1° 

A1/2 

(51) 

/APmin\  _  2 13/2n3/2 

VCvmJ^-  tfA3'2 

(52) 

Comparing  the  twice-minimized  resistance  (52)  with  the  AP  defini¬ 
tion  (32)  that  led  to  the  G  function  (39),  we  conclude  that  the  G 
value  of  the  tree-tree  design  of  Fig.  9  is 


^'3t2'3 


(53) 


In  this  geometric  expression,  the  specified  volume  and  design  of  the 
module  dictate  the  values  of  A  and  t.  For  lower  G  values,  the  area 
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fraction  <j>  should  be  larger,  but  <j>  cannot  increase  indefinitely 
because  it  must  be  smaller  than  1.  Finally,  Eq.  (53)  suggests  that  a 
smaller  n  is  better.  This  brings  us  to  the  question  of  how  to  select 
the  spacing  (L/n)  between  the  D2  channels.  The  answer  is  that  the 
L/n  spacing  must  be  small  enough  (not  greater  than  a  threshold) 
such  that  the  thermal  resistance  to  lateral  heat  conduction  along 
the  solid  plate  (ks),  from  the  D2  channel  to  the  interstice,  is  not 
greater  than  the  thermal  resistance  in  the  same  direction  through 
the  fluid,  from  the  m/n  fluid  to  the  internal  surface  of  the  D2 
channel: 


L/n 

ksD2 


D2 

kJh 


(54) 


In  view  of  the  results  (44),  (50)  and  (51),  the  inequality  (54) 
becomes  a  lower  bound  for  the  projected  tree-tree  area  fraction 

2*|  (55) 


If  tj)c  ~  kf/ks  is  the  order  of  magnitude  of  this  lower  bound,  and  if  </> 
~  </;c,  the  spacing  5  between  D2  channels  should  have  this  order  of 
magnitude: 


By  using  n  ~  L/6  and  Eqs.  (44)  and  (49),  we  find  that  cj>  ~  <j>c  holds 
for  any  n  value. 


6.  Illustrative  example 

As  an  illustration  on  how  to  configure  the  flow  network,  we 
start  with  the  following  assumptions: 

•  The  dimensions  of  each  plate  are  fixed  (Lp  x  Hp),  together  with 
the  area  on  which  the  flow  network  is  embedded  (L  x  ff). 

•  The  cross  section  of  the  channels  is  round.  The  channel  diame¬ 
ter  is  Dj  and  cannot  exceed  Dlimax. 

•  We  learned  from  the  constructal  law  that  vascularization  is  a 
feature  that  allows  us  to  decrease  the  overall  pumping  power 
of  flow  architectures.  Therefore,  channels  with  a  diameter  D2 
will  be  added,  and  their  maximum  value  is  also  Dlmax. 

•  The  total  volume  of  fluid  Vt  is  also  fixed,  which  means  that  the 
plate  porosity  is  fixed  ( cp  =  Vf/V,  where  V  is  the  plate  volume). 

•  The  mass  flow  rate  is  also  constrained. 

Four  configurations  are  studied:  Configuration  1  corresponds  to 
the  sketch  presented  in  Fig.  9.  In  this  case,  the  channels  of  diameter 
D,  serve  for  the  distribution  and  the  collection  of  the  fluid.  Their 
length  L  is  greater  than  H,  which  is  the  length  of  the  n  tubes  of 
diameter  D2.  The  spacing  between  the  channels  of  diameter  D2  is  d. 

Configuration  2  (not  presented  here)  is  also  a  vascularized  net¬ 
work  but  differs  from  configuration  1  in  this  respect:  the  distribu- 


m 


L 


H 


m 


Fig.  10.  Serpentine  covering  the  entire  plate  area. 


tion  and  collection  channels  of  diameter  D,  are  positioned  on  the  H 
size  of  the  domain,  while  the  m  parallel  channels  of  diameter  D2 
are  placed  in  the  L  direction.  These  m  channels  are  separated  by 
a  spacing  t. 

In  Configuration  3,  we  replace  the  vascularized  architecture  by  a 
single  serpentine  of  diameter  Dt.  This  configuration  is  sketched  in 
Fig.  10.  The  spacing  between  the  portions  of  the  channel  aligned 
with  H  remains  at  d. 

Finally,  Configuration  4  (not  shown  here)  is  the  equivalent  of 
configuration  3  but  the  meanders  of  the  serpentine  are  placed 
along  the  L  direction.  In  this  way,  the  spacing  between  the  portions 
of  the  channel  aligned  along  L  remains  equal  to  t,  as  in  Configura¬ 
tion  2.  In  sum,  Configurations  3  and  4  offer  the  opportunity  to  com¬ 
pare  the  merit  of  a  vascularized  network  in  comparison  with  a 
linear  occupation  of  the  area. 

Next,  we  evaluate  the  overall  pressure  losses  in  each  configura¬ 
tion,  assuming  that  the  svelteness  of  the  network  defined  as 
Sv  =  ( HL.yl2IVf 13  is  high  enough  so  that  the  local  pressure  losses 
can  be  neglected.  We  obtain  the  following  results: 

Configuration  1 


For  large  n,  we  have  X)Li  (1  -  £)  =  §  ,  and 


An  ^  .  (n  d  H  \ 

(58) 

V>  =  |(  2D3L+nD22H) 

(59) 

L=d(n  - 1) 

(60) 

Configuration  2 

AP  =  Cvm  — —4') 

\2  D?  mD\ ) 

(61) 

Vf=^(2D]H+mD22L) 

(62) 

H  =  t(m  -  1) 

(63) 

Configuration  3 

AP  =  Cvm ^ 

0f 

(64) 

Vf  =  7f~nH 

(65) 

Configuration  4 

AP  =  Cvm  ^ 

D\ 

(66) 

,/  nDi  , 

Vf  =  ^mL 

(67) 

Recall  that  the  plate  porosity  tp  is  fixed.  In  Configurations  1  and  2, 
we  are  looking  for  the  number  of  channels  to  implement  (n  in  Con¬ 
figuration  1,  m  in  configuration  2)  and  for  the  diameter  ratio  r  =  D,/ 
D2  such  that  the  overall  pressure  losses  are  minimum.  Indeed,  by 
minimizing  the  pressure  losses  with  a  fixed  mass  flow  rate,  means 
to  minimize  the  pumping  power. 

We  assume  that  the  projected  area  of  the  flow  network  has  the 
dimensions  L=140mm  and  H  =  110  mm.  The  plate  porosity  is 
cp  =  0.30.  Water  at  20  °C  flows  at  m  =  9  g/s  through  the  channels 
which  maximum  diameter  is  Dlmax  =  2  mm.  The  results  corre¬ 
sponding  to  the  first  configuration  are  presented  in  Fig.  11.  The 
overall  pressure  losses  are  minimum  when  the  network  is  made 
of  n  =  42  channels  of  diameter  D2  connected  to  the  distribution 
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Fig.  11.  Pressure  losses  as  a  function  of  the  number  of  parallel  channels.  Configuration  1. 


Fig.  12.  Pressure  losses  as  a  function  of  the  number  of  parallel  channels.  Configuration  2. 


channel  of  diameter  Di  and  to  the  collection  channel  of  diameter 
Dj.  This  configuration  corresponds  to  D,  =  D2  =  2  mm,  namely 
r=  1.  We  obtained  AP=  1703  Pa.  Because  the  porosity  is  fixed,  an 
increase  in  the  number  of  channels  corresponds  to  an  increase  in 
r  and  a  decrease  in  the  diameter  D2,  which  leads  therefore  to  an 
increase  in  the  pressure  losses.  Note  that  without  constraining  D, 
and  D2  to  Hi  .max,  the  pressure  losses  would  continuously  increase 
with  the  number  of  channels  in  parallel. 

Fig.  12  shows  the  evolution  of  the  pressures  losses  in  configura¬ 
tion  2.  Again,  a  minimum  can  be  reached  when  m  =  34  channel  in 
parallel.  The  diameter  ratio  is  r=  1.01,  the  diameter  D i  is  2  mm, 
and  the  corresponding  pressure  losses  are  AP  =  1396  Pa.  Note  that 
in  this  configuration  the  number  of  channels  in  parallel  is  smaller, 
yet  the  channels  are  longer,  compared  to  configuration  1. 

Configuration  3,  which  is  the  serpentine  analogous  to  the  vascu¬ 
larized  architecture  of  configuration  1  was  analyzed  in  two  ways. 
First  we  fixed  the  spacing  d  to  be  the  same  as  the  one  obtained 
by  optimizing  configuration  1.  In  this  case,  because  the  porosity 
of  the  plate  is  fixed,  the  diameter  D,  should  be  of  2.92  mm,  which 
is  not  acceptable  because  its  maximum  value  should  be  2  mm. 
Next,  the  diameter  value  is  constrained  to  2  mm.  The  correspond¬ 


ing  spacing  is  d  =  1 .59  mm  for  a  pressure  drop  of  225,01 7  Pa.  In  the 
case  of  configuration  4,  an  equivalent  reasoning  lead  to 
t=  3.23  mm  and  AP=  1.113  x  10s  Pa. 

These  results  confirm  the  superiority  of  the  constructal  design: 
vascularization  is  the  way  to  design  the  flow  network. 


7.  Concluding  remarks 

The  chief  conclusions  of  this  work  are: 

•  The  maximum  power  that  can  be  obtained  from  a  TEM  system 
is  proportional  to  the  temperature  difference  squared. 

•  Counter-intuitively,  parallel  flows  provide  better  performance 
(i.e.  power  density)  than  counter  flows. 

•  Constructal  theory  allows  to  find  the  optimal  geometric  fea¬ 
tures  for  minimum  pumping  power  and  maximum  power 
density 

•  Vascular  architectures  lead  to  better  designs  even  when  the 
maximum  channels  diameters  are  fixed  provided  the  svelteness 
of  the  network  is  high  enough. 
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